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Abstract. In a recent work difference equations (Laguerre-Freud equations) 
for the bi-orthogonal polynomials and related quantities corresponding to the 
weight on the unit circle w{z) = H^iC-^ ~ were derived. Here it is 

shown that in the case m = 3 these difference equations, when applied to 
the calculation of the underlying U{N) average, reduce to a coupled system 
identifiable with that obtained by Adler and van Moerbeke using methods of 
the Toeplitz lattice and Virasoro constraints. Moreover it is shown that this 
coupled system can be reduced to yield the discrete fifth Painleve equation 
dPy a.s it occurs in the theory of the sixth Painleve system. Methods based 
on affine Weyl group symmetries of Backlund transformations have previously 
yielded the dPy equation but with different parameters for the same problem. 
We find the explicit mapping between the two forms. Applications of our 
results are made to give recurrences for the gap probabilities and moments in 
the circular unitary ensemble of random matrices, and to the diagonal spin- 
spin correlation function of the square lattice Ising model. 



1. Introduction 

In a recent work jl2j we undertook a study into differential and difference struc- 
tures associated with bi-orthogonal polynomials for the weight on the unit circle 

m 

(l.I) w{z) = l[{z-z,{t)r^, z = e'^ ee{-7T,7r]. 

i=i 

Our motivation was to eventually use these results to provide recurrences in N for 

the random matrix average 

(1.2) 

1 = 1 '^^^> ^^^> ^"-J-^ J-^ 1=1 l<j<k<N 

Here U{N) denotes the unitary group oi N x N matrices, and the right hand side 
exhibits the corresponding eigenvalue probability density function in the case of the 
Haar (uniform) measure. 

In this work we fulfill our original goal in the case to = 3 by specialising the 
results of ^21 as relevant to the calculation of (|1.2|) with the piecewise continuous 
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version of 1)1. l|l in the case m = 3, 

(1.3) w{z) = t-''z-^-'^{l + zf^' (1 + tzf 



l-i 6* G (tt - tt) 



where /.t, wi, L02 are complex parameters [lo ~ oji -\- 1102) and ^, t = e*"^ are complex 
variables (0 G [0, 27r)). By so studying H1.3() we are also completing a study begun 
in |13j on recurrences satisfied by U{N) averages in random matrix theory from the 
viewpoint of orthogonal polynomial theory. 

The average (|1.2|l with weight (|1.3|) occurs in a variety of problems from math- 
ematical physics. The identification of this average as a r-function for the Painlcve 
VI system 11 has led to its characterisation in terms of a solution of the cr-form of 
the Painleve VI equation. Regarding specific examples, we first mention that the 
case (cji, ^2, y) = (0, 0, 0) is the generating function for the probabihty that the in- 
terval {n — (f), tt) contains exactly k eigenvalues in Dyson's circular unitary ensemble 
CUEat (which is equivalent to the unitary group with Haar measure). Similarly the 
case (a'i,a;2,/i) = (1,0,1) is (apart from a simple factor) the generating function 
for the probability density function of the event that two eigenvalues from CUEAr^2 
are an angle apart with exactly k eigenvalues in between. The case — 2, lu2 — 0, 
fi ^ wi = 1/2 corresponds to the density matrix for the impenetrable Bose gas in 
periodic boundary conditions. It was studied in detail in ^j. Furthermore, in the 
case = 0, one sees that H1.3|) includes as special cases 



JV 



(1.4) ( n ^y'li + + k-'z,y^^) i/k' < 1, 

1=1 ' ' 

JV 

(1-5) (l[{l + l/zir'{l + q'zir) q^<l. 

1=1 ^ ' 

The average H1.4|l is equivalent to the Toeplitz determinant given by Onsager for 
the diagonal spin-spin correlation in the two-dimensional Ising model As a 

PVI T- function it has been studied in and ^Jj. The average H1.5|l occurs as 
a cumulative probability density in the study of processes relating to increasing 
subsequences 

In the case ^ = of (|1.3|l , the same problem as we are addressing here has been 
previously considered by Adler and van Moerbeke , using methods of the Toeplitz 
lattice and Virasoro constraints. The difference equations obtained there were 
not identified with known integrable difference equations. Here we find that our 
formalism of bi-orthogonal polynomials leads to the very same difference equations. 
Moreover, we are able to show that they can be reduced to examples of the discrete 
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Painleve V equation dPy , 

, ifn + 1 - a2){fn + 1 - ao - aa) 



fnifn + as) 



(1-6) gn+ign = t- 

(1-7) fn + fn-1 = . 

gn - 1 gn-t 

where ai > ai + 1, aa i—s- a2 — 1, 0^4 >— > 0^4 + 1 as n >— > n + 1, fundamental in the 
theory of the Pyi system for its relationship to the Backlund transformations ^I] , 
and its association with the degeneration of the rational surface D^P — > o'^p in the 
space of initial conditions PH] , 

The dPy system in relation to U{N) averages (specifically (|1.5(l 'l was first found 
in the work of Borodin |n|,[Z|- Subsequently the present authors 10 used methods 
based on the afhne Weyl group symmetries of Backlund transformations for the 
Pyi system to give dPy recurrences for H1.2(l with weight (|1.3(l . These are different 
to the form of dPy found from our transformation of the recurrences from the 
bi-orthogonal polynomial theory, thus raising the question as to the relationship 
between the two. This we answer by providing the explicit transformation formulae. 

A practical side of our work is that the reduced (in order) form of the Adler and 
van Moerbeke type recurrences are typically much simpler than their equivalent 
form expressed in terms of dPy . We make explicit these recurrences in the cases 
of the characteristic polynomial and gap probability for the CUEjv, and the case 
of the diagonal spin-spin correlation function for the square lattice Ising model. 

In Section 2 general formulae from l2| required in subsequent sections are 
recalled. These formulae are used in Section 3 to derive A^-recurrences in the case 
of weight H1.3|l . and furthermore a transformation to the dPy system (|1.6ll . (|1.7|) is 
given. The explicit transformation between the latter, and the dPy system found 
in relation to as a consequence of the Okamoto r-function theory JHIi is 

established in Section 4. Application of the recurrences to random matrices and 
the Ising model is given in Section 5. 

2. Bl-ORTHOGONAL POLYNOMIAL FORMALISM 

From the viewpoint of our work [12) . the weight 1)1. 3|l is a particular example 
of the regular semi-classical class p.lfl , characterised by a special structure of their 
logarithmic derivatives 

Here V{z), W{z) are polynomials with degV{z) < TO,degT4^(z) — m. We define 
bi-orthogonal polynomials {(/)„(2), (/)„(z)}^o with respect to the weight w{z) on 
the unit circle by the orthogonality relation 

(2.2) /_^^„(C)0„(C)^„(C)=5. 



_ m.ni 
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where T denotes the unit chxle |C| — 1 with C, — e'^, Q G (— 7r,7r]. Notwithstanding 
the notation, is not in general equal to the complex conjugate of 0„. We set 

0„(z) = K„z" + Z„z"-i + m„z"-2 _^ ^^(0)^ 
^„(z) = K„z" + r„z"-i + m„z"-2 + . . . + 0„(O), 

where again Z„, m„, (/>„(0) are not in general equal to the corresponding complex 
conjugate. The coefficients are related by many coupled equations, two of the 
simplest being 



(2.3) 



= +0„(O)0„(O), 



Denote the U{N) average (ll.2|l by In[w]- ^ basic fact that /iv[t(;] can also 
be written as the ToepUtz determinant 



27r< 



0<j,k<n-l 



(2.4) J„H=det 
With the so-called reflection coefficients specified by 
f25) r f 

[■^■•-') ' n — Tin — ) 

it is a well known result in the theory of Toeplitz determinants that 

(2.6) In.Mln M ^ , _ 

Introduce the reciprocal polynomial (t>niz) of the nth degree polynomial (pniz) by 

(2.7) C(z) := z'>„(l/;s). 
Fundamental to our study ^21 is the matrix 

(j)n{z) e„{z)/w{z) \ 



(2.8) 


Yn{z;t) 


where 




(2.9) 


en{z) / 




Jt 


(2.10) 


e;(z) ~z 



h 2niC(^- z 
Thus we obtained the difference system 



(2.11) 



J_ / Hn+lZ (/>«+! (0) 
l^n W«+l(0)z K„+i 
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and the differential system 



(2.12) ~r'Yn • — AnYn 

dz 



1 



W{z) 



i7„(z) + y(z)-^ze„(z) 

Kjri. 



(0) 



e„(z) 



ze*(z) 



ni{z)-V{z)-^Ql{z) 



\ n^n n^n / 

For regular semi-classical weights (jl.lll . the functions 8„(z), 8* (2;), ri„(z) 
and ^n{^) l|2.12(l are polynomials of degree degJl„(z) = degri*(z) = to — 1, 
deg9„(z) = degO* (z) = to — 2, independent of n. Explicitly, for large z 

m 

(2.13) e„(z) = (n + l + yp,)^z"-2 

J = l 

m m m m 3 , 

+ {-K- + i + E^.)E^.-E/'.-.]^ + (- + 2 + E/'.0;^r^|^ 



-(Jl + ^p,) 



k;^ 0«+2(o) 
(0) 

(O)0„(O) ^KnWl l,m-3 ^ 0(2™""*) 



(2.14) r!„(z) = (l + i/2yp,)^'""' 



V2(E ) (E ) + V2 E - E ■ 



i=i i=i 



(n + 2 + y Pj) 



l^n 0ri+2(O) ^"+1 ^ 2™-2 _|_ Q(_2™~3) 



(0) 



Also obtained in |12) , and relevant to the present study, are the small z expansions 

.(0) 



(2.15) e„(z) = [2F(0) - nW^'(O)] 



0n+l(O) 



+ <j [2V'{Q) - i,^M^"(0)]-^^ + [2V{Q) - (n - 1)T4^'(0)] ''"'^""i^^) 



0n+l(O) K„_i(/)„+i(0) 

[{n + l)iy'(O) - 2^(0)]-^ - [{n - l)W^'(O) - 2^(0)]^) /"^"j^, 1 



Kn+1/ '?!'n+l(0) J ^ 

+ 0(z2), 



(2.16) f7„(z) = \/(0) -nVF'(O) 



+ \v'{0) - |7iW^"(0) + (v{Q)^ + [y(0) - nM^'(O)]^ 



+ [1/(0) - nVK'(O)]— - [1^(0) - (n + l)M/'(0)]-^ ^ - ' ^'-^ 



0«+i(O) 
z + 0(z"). 
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Analogous formulae hold for O* (z) and ^n{z) (these can be found in |12j). 

The coefficient functions r2„(z), fi* (z), 0„(z), 6^(z) satisfy sets of difference and 
functional relations given in 12 . Relevant for the present study are the coupled 
equations 

V 0„(O) K„ / 

, K„0„+2(O) K„_l(/)„+l(0) (0) ly (2) _ 
1 TTTT^t^n+U^j , /n^ ^t^ri-ll^^j , ^ ^' 

V l^n (Pn[0) J 
, K„<^„+2(0) , - K„_l(^„+l(0) K„+l W^(z) _ 



Kn+10n+l(O) " Krt(/>n(0) 



(2.19) f7:^,(z) + f^„(z) - + 1^^) e:+i(z) 

\K„+1 (0) / 



(ze„(z) - e„(z)) = 0, 



(2.20) f7„(z) - r2„+i(z) + 1 + z e„+i(z) 

\ Kn+1 Kn+2 / 

+ ^!i±iM^.e„(z) - ^e:(z) = 0, 

•^^-^l) -TTTTrQ^W z0n-i(2:) = — ze„(z) 6„-i(^), 

9)n(0) K„_i (Pn{0) Kn-1 

(2.22) o;(z) - n,,{z) = -^(ze„(z) - e:(z)) + n^^^^ 



Kn Z 



0n+2(O) , , K„+i 

e„+i(z) + e„(z) 



VF(z) 



(2.23) f]:(z) + f7„(z) , rm"'— 

In addition to the above coupled equations, evaluations of the coefficient functions 
at the singular points satisfy bilinear relations jl2| 

K„(/)„+2(0) 



(2.24) 




(2.25) 




(2.26) 





Kn+l(/'n+l(0) 
K„0„+2(O) 



zje„(zj)e„+i(zj) + V (zj), 



l ^n-l 0n+l(O ) 



^n-l{Zj) ^ rn\ ^n{Zj) 



K„+10„+l(O) 

0„+i(O)^„(O)^ , , , , ^ 



^''--■^^"^ ^ "■■e„(z,)e:_i(z,) + F'(z,), 
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(2.27) 



0«+i(O) 

4 MO) 



K„_i(/)„+i(O)0„(O) 2r^*t ^r^ / \ , ^ 



(2.28) ^^^±i#^z,e„(z,)e:(z,) + y^(z,) 



(2.29) 



It is these relations that lead directly to one of the pair of coupled discrete Painleve 
equations. Deformation derivatives of the system H2.8() with respect to arbitrary 
trajectories of the singularities Zj{t) are given in [T2|, leading to the Schlesinger 
equations from the theory of isomonodromic deformations of linear systems. The 
particular result that we require here is 



(2.30) 
(2.31) 



r„ dt 



1 d 



3 = 1 



zjdt^' V{zj) 
1 d nt,_^{z,) + V{z,) 



3 = 1 



dt 



V{z,) 



3. The Semi-classical Class m = 3 and dPy 

3.1. Coupled Recurrences for the Reflection Coeflticients. Here we will con- 
sider the general theory from '12' revised above, applied to the simplest instance 
of the semi-classical weight Hl.l|l . namely m = 3 singular points with two fixed at 
z = 0,-1 and the third variable at z = —1/t. Explicitly we consider the unitary 
group average (|1.2|) in the case of the weight p.3|) . For this to make immediate 
sense we require t G T, but it can analytically continued off the unit circle. Note 
that when t G T, ^,uji,uj2, £, G K and ^ < 1 the weight (|1.3|l is real and positive. 
The corresponding Toeplitz determinant (|2.4(l is then hermitian and as a conse- 
quence f„ is the complex conjugate of r„, but generally this is not the case. For 
parameters such that K(/i), 5ft(ti'i) > —1/2, t <E T the Toeplitz matrix elements, 
Wj-k say, can be evaluated in terms of the Gauss hypergeometric function and so 
analytically continued into the general parameter space. There are several forms 
for this, which exhibit the analytic structure about the special points t = 0, 1, oo. 
We will make note of two such expansions, relating to the special points t — and 
t = 1. 

Lemma 3.1. For all points in parameter space that the functions below have mean- 
ing, the analytic continuation of the Toeplitz matrix elements w„ for the weight U.<-i\) 
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is given by 
(3.1) 



r(2cJi + 1) 



(— 2^, —n - ^ — L0] \ - n - ^ + uj]t) 



r(l + n + /i + w)r(l — n — /i + ti-") 

g ±M^-^) + i)r(2^i + 1) 

27ri r(2/i + 2cji + 2) 

X 2Fi(2^ + l,l + n + ^ + w;2/i + 2wi + 2;1 -t), 

where the ± si^n is taken accordingly as Ini(i) ^ 0. This can also be written as 

(3.2) rw„==n+^— — — ^ ^ r "FTT^ — ; — \ — vfti 

2z sin7r(n + /i — cj) J L (1 + n + + uj)l [1 — n — fi + uj) 
X 2^i(~2/i, — ri — /i — cj; 1 — ri — /i + tj; t) 

g±,r»(n+M-<S) r(2/i+l) 

2i sin7r(n + jj, — Co) T{1 + n + /i — tL')r(l — ?i + /i + tj) 



X t 



(1-0 



2/^+2^1 + 1 



i^i(2^ + 1, 1 + n + /i + tj; 1 + n + ^ - cj; t). 



Proof. This follows from the generalisation of the Euler integral for the Gauss 
hypergeometric function and consideration of the consistent phases for the branch 
cuts linking the singular points, see PP- 91i section 17 " Verallgemeinerung der 
Eulersche Integrale" . □ 

Remark 3.1. The first of these forms was given in jl()| . 

For the description of (|1.3|l in terms of the characterisation 1)2. l|l of general 
semi-classical weights we have 

(3.3) m = 3, = {0,-1,-1/0, {p,}%, = {-ii-u,2u;,,2^i}, 

1-t 



(3.4) 



2V{0) = -{ji + Lj)t-\ W'{0) = t-\ V{-t-^) = n- 



Also, according to H2.13)l - I|2.16|l . the explicit forms of the coefficient functions in 
this case is 



(3.5) Qn{z) 

(3.6) e^(z) = 



UN 

kn+1 



{N + 1 + fi + u})z 



rN+i 



{N + fi + Lu)t-^ 



■TN+l 



(3.7) nN{z) = [l + ^{fi + Ld)]z^ 



+ \ (iV + 2 + ^ + tZ.)(l -rAT+ifjv+i) 



rN+2 In+i 



rN+i KN+i 



+ [l + i(/i + C^)] 



1 + t 



z~ i(^ + w)]i~\ 
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(3.8) n*j,{z) = -^{ii + Q)z'' 

+ i [N + 11 + Lo)[l ~ TN+irM+i)- ^{^l + uJ)— Vuji + - > z 

L kat+i »"a+i t t } 

We are seeking a closed system of recurrences for r^r and fjv, since according 
to (|2.6(l these quantities determine the Toephtz determinant, or equivalently the 
U{N) average. One such recurrence is quite straight forward. 

Lemma 3.2. The reflection coefficients for the weight satisfy the homoge- 

neous second- order difference equation 

(3.9) (iV + 1 + ^ + u!)trN+ifN - (iV - 1 + ^ + (I>)irArfAr_i 

= {N + 1 + fi + u;)fN+irN - {N — 1 + ^ + uo)fMrN-i- 

Proof. This equation can be found immediately from the general theory of Section 
2 in many ways. By equating coefficients of z in the functional-difference equation 
(|2.19|l using H3.5I I3.6I3.7I3.8|I . all are trivially satisfied except for the z coefficient, 
which is precisely 1)3. 9|l . Similarly starting with H2.2()|l and employing 1)3.51 13.6I3.8|I . 
one finds 1)3.9(1 . Alternatively one could start with either ()2.21)l or ()2.23(l and arrive 
at the same result □ 

We will use this result in the derivation of a sequence of lemmas, which lead us 
to the sought closed system of coupled recurrences. 

Corollary 3.1. The sub-leading coefficients In, In satisfy the linear inhomogeneous 
equation 

(3.10) (iV + /I + Lu)tlN -{N + + uj)In = N [^{t - 1) + w - wt] kn- 

Proof. By substituting the general expression for the first difference of In, In using 
()2.3(l in 1)3.9)1 one finds that it can be summed exactly to yield 

(3.11) (7V+l + ^ + w)i-^-(7V + l + M + t^)-^ 

KN+l KN+1 

- {N + uj)t— + {N + fi + Lu)— = fi{t - 1) + u) - ut. 

KN KN 

This can be summed once more to yield the stated result. □ 

Remark 3.2. One could alternatively proceed via the Freud approach (see also 
and consider the integral 



(3.12) / + + 



/T 



fi + Lo 2loi 2fit 



l + z l + tz 



■w{z)(j)N{z)(i)N{z). 
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Here we recognise the logarithmic derivative of the weight function in the integrand 

(3.13) w^^_f^^2uJ^^J^^ 
^ ^ w z 1 + z l + tz 

and by evaluating the integral in the two ways we find a linear equation for In, 
namely (|3.11() . 

Lemma 3.3. The sub-leading coefficients are related to the reflection coefficients 
by 

In/ UN + tlN/nN - N{t + 1) 

(3.14) = ^ " '"^''^ [(AT + 1 + ^ + Oj)trN+i + {N-l + n + uj)rN-i] , 

rN 

(3.15) = ^ ~ '^^''^ [(AT + 1 + ^ + uj)fN+i + (iV - 1 + M + i:')trN-^i] ■ 

rN 

Proof. The first relation follows from a comparison of the coefficients of z for flN{z) 
given the two distinct expansions, the first by (|2.14|) which reduces to (|3.7|l and the 
second by the specialisation of H2.16|l . The second relation follows from identical 
arguments applied to f2^(z) or by employing H3.9|l in the first relation. □ 

Remark 3.3. The first result appears in the Magnus derivation (2(1) for the gener- 
alised Jacobi weight, with 9i = tt — 0, 02 — t^, ol = /x, /3 = 0^1,7 = —loi. Then 
Equation (14) of that work is precisely (|3.14l) . 

Remark 3.4. The Magnus relation (|3.14() can also be found by employing the Freud 
method. In this one uses integration by parts on the integral 

f dz 

(3.16) / ——z~\l + z){l + tz)w'{z)(bN+i{z)(l)N{z), 

and in the term involving 0'^^-^(z) one employs (I2.12|) for the derivative and H3.5|l . 
()3.7|) for the coefficient functions. Equating this expression to a direct evaluation 
of the integral then yields (|3.14() . 

Lemma 3.4. The sub-leading coefficient In can be expressed in terms of the reflec- 
tion coefficients in the following ways 

(3.17) 2t^ =(iV + 1 + M + cS)<(^ - rN+ifN) + (iV - 1 + /i + cj)^ 

KN rN rN 

-{N -1 + fi-h <Z>)rjvfjv_i -\- {N -\- fi - uj)t + N - fi -\- Lu, 

(3.18) ={N + 1 + M + Lo)^^^ + {N-l + ^I + Ld)tC-^ - rjvfAT-i) 

rN rN 

~{N + 1 + 11 + uj)trN+ifN + {N + fi - uj)t + N - ^ + Q, 
as well as analogous expressions for In- 
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Proof. The first expression follows from a comparison of the z'^ coefficients for 
Qn{z) evaluated using both H2.13|l and l|2.15|l . The second relation follows from an 
applying the same reasoning to Q*j^(z). □ 

We will refer to the order of a system of coupled difference equations with two 
variables, r„,f„ say, as q/p where q S Z>o refers to the order of r„ and p G Z>o 
refers to the order of f„ . 

Corollary 3.2. The reflection coefficients of the OPS for the weight iLf-^) satisfy 
the 2/2 order recurrence relations 

trNTN-l + TN-lfN - t - 1 

1 - rAffjv 



[(iV + 1 + /i + u))trN+i + (TV - 1 + /I + uj)rN- 



(3.19) -- — ^ [(A^ + M + w)fjv + (iV - 2 + /i + uj)trN-2\ , 

= ^ [{N + 1 + n + uj)fN+i + {N -1 + n + id)tfN-i] 

rN 

(3.20) ^ l-rN-irN-i [(^ ^ ^ 4, ^)^^^ + {N - 2 + fi + Lo)rN-2] , 

rN-i 

and those specifying the solution for 1^2. b]) have the initial values 

(3.21) To = To = 1, ri = -w_i/wo, ri = -wi/wo, 
where the Toeplitz matrix elements are given in yS. 2\) . 

Proof. Solving (|3.14|) for the combination oUnjIn and differencing this, one arrives 
at (|3.19() . This however is of order 3/1 but by employing 13.9|l we can reduce the 
order in rjv of the recurrence to second order. The other member of the pair (|3.2U|) 
is found in an identical manner starting with l|3.15|l . □ 

Remark 3.5. The second-order difference ()3.19|) also follows immediately from equat- 
ing the polynomials in z arising in the functional-difference (|2.17l) . after employ- 
ing (|3.5I3.7|) . The other member of the pair, (|3.2U|) . follows from the functional- 
difference (|2.18ll . after using (|3.()I3.8|I . 

3.2. Relationship to recurrences of Adler and van Moerbeke. In their most 
general example, Adler and van Moerbeke 1 also considered recurrences from 
the U{N) average with weight H1.3|l (^ = case). Their study proceeded via 
the viewpoint of the Toeplitz lattice and Virasoro constraints. In terms of their 
variables we should set Pi ~ P2 = 0, di = t~^/^,d2 — t^^^, and without loss of 
generality 7" = 73 = 0. For the other parameters j = fj, — lu,j[ = 2ci;i,72 = 2/i. 
There is a slight difference in the dependent variables due to the additional factor of 
t, so that we have the identification xn = (— l)^i^/^rjv, yN — (— l)^t^^/^fAr and 
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WAT = 1 — rNTN- Generalising their working one finds that their Equation (0.0.14) 
impUes 

(3.22) - (TV + 1 + /i + u))xN+iyN + {N + l + ^l + uj)xNyN+i 

+ (iV - 1 + /i + u})xNyN-i - (A^ - 1 + /i + uj)xN-iyN — 0. 

Now by transforming to our r^v , tn and employing 1)2. 3|l one finds this is precisely 
(|3.9|l . which we showed is solved by (|3.10|) . Their inhomogeneous Equation (0.0.15) 
now takes the form 

(3.23) -vn[{N + 1 + fi + uj)xN+iyN-i + N + ^ + uj] 

+ vn-1 [{N ~2 + fi + ui)xNyN-2 + N-l + fi + uj] 

+ XNyN-lixNVN-l + t^^^ + t^^^^) 
= -VI [(2 + ^ + U))X2 + 1+ fJ, + Lu]+Xl{xi+ t^/^ + t"^/^). 

Upon recasting this into our variables and manipulating, it then becomes 

trwrw^i + rNrN-i - t-1 

- ^ ~ ^^^'^ [{N + 1 + n + Lj)trN+i + (N -1 + 1^1 + uj)rN-i] 
+ 1 - ^^-i^^-i [(AT + ^ + u:)rN + (iV - 2 + M + cD)ifjv-2] 

_ 1 - (1 - rifi) [(2 + /i + U))tr2 + 1 + + ^] + ri {tn - t - I) 

However using the identity 
C2fi(a, c; a;) 

= \c+{l + b-a)x\2Fi{a,b+l]c+l;x) - -^-^(1 + c - a)a:2i^i(a, 5 + 2; c + 2; a;), 

c + 1 

we note that the right-hand side is identically zero for the initial conditions (|3.21|l 
and the recurrence is not genuinely inhomogeneous, thus yielding our first relation 
above, (|TTn|) . 

3.3. Transformation of the recurrences to dPy • We seek recurrences for 
rjsijfiy which are of the form of the discrete Painleve system (|1.6|l . p. 7(1 . For 
this purpose a number of distinct forms of the former will be presented. 
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Proposition 3.1. The reflection coefficients satisfy a system of a 2/0 order recur- 
rence relation 

(3.24) |(1 - rNfN) [{N + l + fi + oj){N + n + oj)trN+i 

-{N + /i + w)(iV - 1 + /i + uj)rN^i] + N{N + 2uJi){t - l)rN^ 
X |(1 - rNfN) [{N + l + + Q){N + fi + u)trN+i 
-{N + + uj){N - 1 + /i + u})rN-i] + {N + 2/i)(7V + 2/i + 2uji){t - l)rN^ 
= -{2N + 2/i + 2wi)2i(l - rNfN) 
x[{N + 1+ fi + Lj)rN+i + {N + ^ + Lu)rN] 

x[{N + ^ + uj)rN + (iV - 1 + ^ + uj)rN-i] , 

and a 0/2 order recurrence relation which is just \S.24^ with the replacements lo uj 
and t^^/'^r^ ^ t=Fi/2^^. 

Proof. Consider first the specialisation of (|2.24|) to our weight at hand at the sin- 
gular point z = — 1, and we have 



(3.25) I - Nt-' - (TV + 1 + M + c^)^^ + - t-^)\ 

iKN Km rN J 



+ 



"N 



N + fl + Ld + 



N rN 
{N-l + fi + Lu) rN-i 



t 



{N + fl + Lu) 



+ (iV + 1 + /i + (I)) 



rN 
rN+i 



rN 



t - 1 



by using (|3.5I3.7|) . Similarly 12.24|l evaluated at z = — 1/i yields 



(3.26) lhL^N-iN + l + f, + u,)^'-^ + ^,{t-' 1)1 
I Km ny rN J 



•^N-l ,-1 



K 



N 



N + iJ. + Q+{N -1 + fi + uj) 



rN-i 



N + ^i + uj+{N + l + fi + u) 



rN 
rN+i 



TN 



t - V 



The first relation follow by eliminating between (|3.25f) and (|3.26|l . whereas the 
second follows from an identical analysis to that employed in the proof of Proposi- 
tion]^] but starting with the bilinear identity l|2.25|) . □ 
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Proposition 3.2. The reflection coefficients also satisfy a system of 1/1 order 
recurrence relations the first of which is 

(3.27) 

I - (A^ + + w)(l - rNfN)t [{N + l + fi + u)rN+ifN + (N - 1 + n + u)rNfN-i] 
+ 2{N + fi + uj)'^rj^f% - {N + fi + Luf{t + l)rNfN ~ 2(A^ + /i + - l)rjvfjv 

X 

I - (TV + ^ + ijj){l - rArfjv)^ [(A^ + 1 + /i + w)rAr+ifAr + (A^ - 1 + ^ + (I')rArfjv_i] 

+ 2(iV + ^ + t^)V^f^ - (TV + /i + ujf{t + l)rArfAr + 2(iV + /.t + w)cj(i - l)rjvfw 

= — [2(A^ + ^ + Lj)rNfN +L0 ~ ujf (1 — r^rN) 
x[{N + l + (1 + uj)trN+i + {N + fi + uj)rN] 

X[{N + H + L0)fN + {N -1 + 11 + U})tfN-l] , 

and the second is obtained from \^.2T{ l with the replacements uj ^ u) and t^^^'^r 
t^^/^f 



3 

3 ■ 



Proof. The specialisation of (|2.26() to the weight (|1.3|l evaluated at the singular 
point z — —1 is 

(3.28) {^L _ ATt-i + [N + ^i + c.)i-i%i + ..1(1 - i-i)|' 

+ [(iV + 1 + + w)rAr+i + (iV + /i + w)t"VAr] 



X [(iV - 1 + /i + (Ii)fjv_i + (iV + ^ + w)t"^fAr] = CjJ 

by using (|3. 513. 7(1 . Similarly H2.26|l evaluated at z = — 1/t yields 
(3.29) {hL^N + {N + ^l + uJ)'^^+^i{t-^-l)y 

^2 



t-1-^ 



+ ^^^-^ [{N +l + ^i + w)rjv+i + (iV + /i + w)rjv 

K 



t-V' 



N_-l 
2 

x[{N -1 + fi + u))fN-i + iN + fi + uj)fN] = 

Again eliminating Ijy between these two equations yields the recurrence relation 
(|3.27|l . The second follows in the same way starting with (|2.27() . □ 
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Proposition 3.3. The reflection coefficients satisfy an alternative system of 1/1 
order recurrence relations the first of which is 

(3.30) {N +1 + fi + Ld){N + fi + Ld)trN+irN 

-{N+1 + fi + uj){N + + uj)fN+irN + (cj - ^i){uj + fi){t - 1) 

X {N + 1 + fi + uj){N + fi + uj)trN+ifN 
- (TV + 1 + /i + uj){N + + uj)fN+irN + {uj- ^i){uJ + ^J){t - 1) 
= {Cj - uj f[{N +1 + + u))trN+i + {N + ^1 + ijj)rN] 

X [(TV + 1 + ^ + L0)fN+l + {N + ^l + Ld)tfN] , 

and the second is again obtained from \3.8U\) with the replacements a; <-> (D and 

Proof. The specialisation of (|2.28() to the weight 1)1. 3|l evaluated at the singular 
point z ~ —\ is 



(3.31) ( + (iV + /i + w)r7v+irjv + + (^ - iuj2)t\ 

= [{N + 1 + 11 + iJj)trN+i + {N + 11 + Loy^] 

X [(TV + 1 + /i + Lo)fN+i + {N + 11 + Ld)tfi^] + ujl{t - if, 



by using (|3. 513. 7(1 . Similarly 12.28|l evaluated at 



-1/t yields 



(3.32) (i^ + (7V + /i + tj)fAr+irAr+(I. + ^A^ 
L KAT+l J 

= t[{N +1 + ^l + Cj)rN+i + {N + ^1 + uj)rN] 

X [(A^ + 1 + ^ + L0)fN+l + {N + ^l + Lu)fN] + fi'^it - 1)^ 

Again eliminating /jv+i between these two equations yields the recurrence relation 
(|3.30|l . The second follows in the same way starting with (|2.29() . □ 

Remark 3.6. Note that the recurrence system (|3.24ll and its partner is quadratic in 
T'N+i, and fAr+i, tat-i, the system (|3.27f) and its partner is also quadratic in 

rN+i,rN-i and fN+i, r^-i, and likewise H3.30|l is quadratic in rjv+i, r^v+i- This 
renders them less useful in practical iterations than the higher order systems that 
are linear in the highest difference. By raising the order of one of the variables by 
one we can obtain a recurrence linear in the highest difference. 
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Corollary 3.3. The reflection coefficients satisfy a system o/a 2/1 order recurrence 
relation 

(3.33) {N + H + - uj)t{l - rNfN)rN+i 

+ {N -l + fi + uj)[2{N + ^ + uj)rNfN Lo]rN-i 
- {N ~ I + fi + uj){2N + 2/1 + 2LJi)tr%fN-i 
+ [{iD-uj)N{t + 1) - (2^ + 2wi)[/i(l - t) + ut - u]]rN = 0, 

and a 1/2 order recurrence relation which is again obtained from 1^3. 3S\) with the 
replacements w a) and t^^^'^rj <-> t^^^'^fj. 

Proof. The solutions for the sub- leading coefficient In, In that arise from the si- 
multaneous solution of H3.28I3.29|) and H3.31I3.32|) respectively are given by 



(3.34) t— ^ \{N + ^i + uj)t{l~rNfN) 

Kn 



h 

Kn 

X [{N + 1 + fi + iD)rN+ifN + {N ~1 + fi + oj)rNfN^i] 
+ {N + ^i + uj)[N{t+ 1) - - t) - ujt + u}]rNfN + {uJ + fi)[fi{l - t)+ujt-uj]^ 

[2(A^ + ^ + u})rNfN + Cj ~ uj]. 



(3.35) 
t 

KN 



^{N + fx + uj)[{N - 1 + ^1 + Lu)trNfN~i -{N -1 + ^1 + uj)fNrN-i] 
+ (lu + /i)[/i(l — t) + Lut — uj]^ ^ {Co — uj) , if (D 7^ 

and the corresponding expression for In/^n under the above replacements. Equat- 
ing these two forms then leads to H3.33|l . □ 

The systems of recurrences that we have found are in fact equivalent to the 
discrete Painleve equation associated with the degeneration of the rational surface 
D^^i^W and we give our first demonstration of this fact here. 

Proposition 3.4. The N -recurrence for the reflection coefficients of the orthogonal 
polynomial system with the weight jJ.^] is governed by either of two systems of 
coupled first order discrete Painleve equations 1^1. b]l . jJ.7| ). This first is 

XfN + N)ifN + N + 2fi) 



(3.36) gN+i9N = t- 



(3.37) fN + fN-i - 2lo, 



/w(/w-2t^i) 

N -1 + fi + u! [N + ^i + uj)t 



gN - 1 9n -t 

subject to the initial conditions 

(3.38) + + ^ + /o = 0. 
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The transformations relating these variables to the reflection coefficients are given 
by 



N - 1 + 1.1 + uj + {N + ^ + u;)- 



(3.39) 



(3.40) 



gN = t- 



rN-1 



N-l + n + uj+{N + fi + Lu)t 



rN-1 



fN = 



1 



1 - 1 

The second system is 



KN rN 



(3.41) 



(3.42) 



Qn+iQn — t 



In + Jn-i — 2^ 



^i ifN + N)ifN + N + 2coi) 

fN{fN~2fi) 

N + fi + LU {N - 1 + fi + uj)t-^ 



gN - I gN - t ^ 

subject to the initial conditions 

^ + (I'+(l+/i + w)f~^fi 



(3.43) 



5i 



/o = 0. 



^ + [I'+(l+^ + Uj)fl 

The transformations relating these variables to the reflection coefficients are given 
by 

TN 



A^-l + ^ + (Ii+(iV + /i + ijjfl-^- 



(3.44) 



(3.45) 



gN 



TN-l 



iV-l + ^ + w+(iV + /i + cj) 



TN 



In = 



1 



1 - t 



-t^IL +m+{N-l + fi + - rNfN)t^^^ 
KN rN 



Proof. Consolidating each of (|3.25|l and (|3.2t)|) into two terms and taking their ratio 
then leads to (|3.36|) after utilising the definitions H3.39I3.40() . The second member 
of the recurrence system (|3.37|) follows from the relation 



(3.46) 



In+1 

HN+l 



In 



= {N+2+^i+uj){l^rN+ifN+i)^^^ + {N+fi+u;)t-^— {N+l+fi+ui)rN+ifN 

rN+i rN+i 

- 2uoi - 2fit^^ + (A^ + 1 + /i + w)(l + i"^), 



rN 



which results from a combination of H3.17|l and H2.3|l . and the definition (|3.4U|) . All 
the results for the second system follow by applying identical reasoning starting 
with lE23l. □ 

3.4. Evaluations in terms of generalised hypergeometric functions. In the 

special case ^ = of the U{N) average H1.2(l with weight H1.3(l . we know from 

that an evaluation in terms of a generalised hypergeometric function 2-^1^"^ 
possible. Let us first recall the definition of the latter. Given a partition k = 
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(ki,K2, . . . ,kn) such that ki > K2 > • • • > kn > one defines the generahsed, 
multi- variable hypergeometric function through a series representation |24l llftj 
(3.47) 

°° [ai]^^ ■ ■ ■[ap]K^ Sf,{ti, . . . ,tN) 



pFq^Hai, ■ ■ ■,ap;bi, 



K>0 



[bi]i'^---hi'^ 



for p,q ^ Z>o- Here the generahsed Pochhammcr symbols are 



N 



(3.48) 

the hook length is 
(3.49) 



j=i 



where a{i, j),l{i, j) are the arm and leg lengths of the (i,j)th box in the Young 
diagram of the partition k, and s^iti,. . . jt^) is the Schur symmetric polynomial 
of N variables. The superscript (1) distinguishes these functions from the single 
variable N — 1 functions and also indicates that they are a special case of a more 
general function parameterised by an arbitrary complex number d ^ 1. 
With this definition recalled, the result of JT] (see also [5|) reads 



N 



(3.50) ( H zr''""(l + ^0'"' (1 + t^if'') 



1=1 



U{N) 



N-l 

" 5 W+^^ + ^ + J)^a-^^ + ^ + J) 

^ 2^1^'' (~2/x, ~n-uj;N - ^ + u;;ti,..., ijv)|ii=...=t„=t, 

subject to 5R(a;i) > — 1/2 and \t\ < 1. The new observation we make here is that 
the reflection coefiicients determining (|3.5U|I can similarly be written in terms of 
the 2-^1^'' function. 

Proposition 3.5. With Ipf[w] in 1^2. b]) given by l^3.5(J\) . the corresponding reflection 
coefficients are given by 

jY (/i + Cj)Ar 



(3.51) rN = {-lY 



(1 - /i + uj)n 



(3.52) fN^i-l) 



„ (-/i + w)Ar 



(1 + ^ + u;)n 

f[^\-2^i, -1- h-lu;N -1- ^i + u;ti,...,tN) 



2F[^^ (-2^, -^i~u\N ~ ^i + uj-ti 



, tN) 



tl . . . 1 
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Proof. For e = 0, ±1 we define the Toeplitz determinants or U{N) averages 



(3.53) /,tH:=det 



From the Szcgo theory we know 



The result now follows from H3.5()|l . □ 
According to the definition (|3.47() . 2^1^'' is normalised so that 

(3.55) 2Pi^Haib;c;ti, . . . ,tN)\ti=...=t,^=a = 1- 

At the special point when each ti equals unity, the analog of the Gauss summation 
gives the gamma function evaluation j24| 

(3.56) ^^i^H^'^l^i -iJ.-uj;N - fi + uj;ti,..., tN)\ti=...=tM=i 

when ^{fi + LUi) > -I/2, 3?(-/i + uj) > -1. 
The result H3.55(l tells us that at t = 

/or-7^ , _/ (M + ^)jV - _( (-^ + '^)jV 

(3.57) rN - (-1) J- T^^' - (-1) 7— ■ 

while (13.561) tells us that at i = 1 



(3.58) 7Ar-(-l) — — TT— , rjv - (-1) — — ■ — . 

(1 + ^ + U!)n (l + fi + UJjN 

We note that the corresponding values of In can be calculated as 
(3 59) jjv ^ __(m+^)A^ In _ it^ + io)N 

KN (N — fl + U))' KN (N + fl + Lo) 

3.5. Two remarks. We conclude this section with two remarks. The first relates 
to the limit transition from the average (|1.2|l with weight H1.3|l to the average 

N 

(3.60) (^l[zl^-^^/'\l + zir+''e''") ^, 

1=1 

studied from the viewpoint of bi-orthogonal polynomials in 13 . The latter can be 
obtained as a degeneration of the former by the replacements a; + /i 1— > z^, w — /i >— > 
fi,t t/2fi and then taking the limit — > 00. The coefficients of the orthogonal 
polynomials tat, / at remain of 0(1) in this limit. Then we see the explicit degenera- 
tion of the following equations - H3.1()|l Equation (4.23) ^HIi the recurrence rela- 
tions Ipn7|l ^ Equation (4.60)1121, -> Equation (4.9) [El and its conjugate 
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to Equation (4.10)IT51 modulo the identity Equation (4.5) 13 , and the hypergeo- 
metric functions H3.50(l Equation (4.24)^, l|3.51|l — > Equation (4.26)^, and 
^ Equation (4.27)[T3]. 
The second remark relates to the choices of parameters, noted below (|1.3|l . for 
which the weight is real and positive, and consequently Tn = Let us suppose 
furthermore that 0^2 = 0, ^ = 0. Then with t = e'"^, 1)1. 3|l reads 

(3.61) w{e''^) = |2cosi6ip"'|2cosi(6' + 0)p'^. 

It follows from this that w^i = twi, and this in H3.21|) tells us that fi = tri. This 
initial value, together with the initial value given by the first equation in H3.21|l . 
allows a structural formula for r„ in this case to be obtained. 

Corollary 3.4. Let U.,'^) be specialised to L'i.ffl]) . Then the reflection coefficient 
Tn, which is related to fn by complex conjugation, has the form rn = f^^'^Xn, where 
the Xn are real. 



Proof. Setting C1J2 = in 1)3.9(1 we note this can be rearranged as 



(3.62) (71 + 1 + /i + uj) 



+ (n - 1 + /i + cj) 



^ '11. 



= 0. 



It is easy to verify that this has the solution f„ = i"r„ which is furthermore 
consistent with the initial conditions. The result now follows from the fact that f„ 
is the complex conjugate of r„. □ 

4. The t-function Theory for Pyi 

In a previous study 10 the average (|1.2() with weight (|1.3|) has been shown to 
satisfy recurrences involving dPy , distinct from those isolated in Proposition 13.41 
Here we will present the explicit transformation between the two systems. We will 
also show how the relationship between the U{N) average and an average over the 
Jacobi unitary ensemble, studied from the viewpoint of r-function theory in [llj . 
can be used to deduce a further characterisation involving dPy . 

The study ^Hl is based on the Okamoto r-function theory of Painleve systems. 
In the case of Pyi , a Hamiltonian H is defined in terms of coordinate and momenta 
variables a time variable t, and parameters ao,ai,a2,a3,ai G C subject to 
the constraint 

(4.1) ao + ai + 2a2 + as + a4 — 1, 
according to 

(4.2) K:=t{t-l)H 

= 9(9 - 1)(9 - t)P^ - [ai{q -l){q-t)+ a3q{q - t) + {ao - l)q{q - l)]p 

+ a2{ai + a2){q - t). 
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The time evolution of q and p is governed by the Hamilton equations 

dq dH dp dH 
^ ■ ^ dt^ 9^' dt " " a^' 

with the relationship to Pyi coming from the fact that eliminating p gives the sixth 
Painleve equation in q. A crucial quantity in the development of this viewpoint 
given in (22) is the r-function, defined so that 

(4.4) H=j^\ogT. 

It was shown in ^I] that the U{N) average (|1.2() with weight 1)1. 3(1 is a r-function, 
for the Pvi system with the parameters 

(4.5) (ao, ai, "2, "3, a^) ^ (^N + 1 + 2uji,N + 2^, -N, - - , 
and so we can write 



(4.6) T^i[iV](t;/x,c^i,c^2;e) 



N 



The algebraic approach used in |10| makes use of a particular shift operator, 
or Schlesinger transformation L, constructed from compositions of fundamental 
reflection operators and Dynkin diagram automorphisms of the Pyi affine Weyl 
symmetry group Wa{D'^'^) = (sq, si, S2, S3, S4, ri, ra, r4). Application of the oper- 
ator L allows for t^^[N] to be computed by a recurrence scheme in N involving 
auxiliary quantities which satisfy the dPy recurrence. 

Explicitly, the operator with the property of incrementing while leaving the 
other parameters unchanged is Lq^ ~ ?'iSoSiS2S3S4S2- On the root system param- 
eters aj in H4.2I) it has the action 

(4.7) : ao ^ ao + l,ai ^ ai + 1, a2 ^ a2 - 1, 

while a3 and a4 remain unchanged. Studying the action of Lqi on r^^ [A^] in the 
context of the Okamoto theory, the following result was obtained in 1^ . 

Proposition 4.1 ( 10 ). Let {gN, .fN}N=o,i...., satisfy the discrete Painleve coupled 
'erence equations associated with the degeneration of the rational surface d'^^^ 



t ifN + N + l)ifN + N+l + ^l + Q) 

(4-s) gN+WN = - — ^-77 ^ , 

(4.9) /jv + /w-i-M + ^ + T + Ti T^^' 

gN-^ t\gN-i-)~gN 

where t — 1/(1 — e"^) subject to the initial conditions 

go = fo = il + ^i + Co){qo + + oj)qo - (2c.i + 1)?°^'^° 



qo-l qo-t 
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with 



(4.10) '^" = Kl%^l°^^''^'^^(^^') 

Define {qN,PN}N=o,i,... by 



9N 



5Af - 1' 



{9N 1? , 

PN = JN 

9N 

-{N + l + fi + Q)^^^^ - (/i + uj){gM - 1) + (A^ + 1 + 2lui)- " ^ 



t+(l-<)(/Ar' 

T/ien ro(e"^) = 1 and Ti{e"^) = wo(e*"^) as given by jiilAij)) . {TV}Ar=2,3,... 
is specified by the recurrence 

(4.11) -(7V + ^ + ^)(7V + ^ + ^)Zim|i^ 

JV 

= gAf(gAf - ^)p% + {2n + 2uJi)qNPN - {p + ^)pn - N{N + 2^ + 2uji). 

An immediate question is the relationship between the Hamiltonian variables 
Qn,Pn in Proposition 14.11 and the reflection coefficients rpf^f^ relating to t^^[N] 
as studied in Section 3. In fact the quantities are linked by systems of equations 
given in the following result. 

Proposition 4.2. The transformations linking the Hamiltonian variables qntPn 
in Proposition \4. 1\ to the reflection coefficients r^Vj^Af for the system of orthogonal 
polynomials with the weight H.f-i]) are given implicitly by 

{N + H + Lu)r]yfN 1 
QNPN + fl + LJ = J— — T = ■ - 

(4.12) X 



(N + 2uji){qN -l)-t— + Nt+{N + l + fi + uj){l - rNfN)t'''^^^ 



(4.13) 



(4.14) 



KN rN 
^ [N + + Ld)[{N + ^i + uj)rNfN - fi + d] 

X - ^Ji . , 

{N + 2uJi)qN +t—-Nt-{N-l + fi + - rArfA,)i^^ 

KN TN 

{qN — 1)pa + M + ^^ = {N + h + uj) [{N + fi + Lj)rNfN — fi + uj] 

gjV - 1 

{N + 2L0i){qN -l)-t— + Nt+{N + l + ^l + u)){l - rNfN)t^-^ 

KN rN 

{N + fi + Lu)rNfN 1 



{N + ^ + uj)rNrN - fi + uj qN 



(4.15) 



{N + 2uji)qN +t— - Nt- {N -1 + fi + uj){l ~ rNfN)t^'^ ^ 



KN rN 
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Proof. We require in addition to the primary shift operator L^i generating the 
1-^- iV + 1 sequence another operator which has the action iLU2 i— > iuj2 — 1 • This is 
the secondary shift operator Tg^^ = riS4S2SoSiS2S4 and has the action Tg^^ : — > 
as + 1, q;4 — > q;4 — 1. From Table 1 of we compute the actions of Tg^^, T34 on 
the Hamihonian to be 

9n - 1 



+ {ao + a4- l)(g'n - 1) - (0^2 + a3)(ai + ^2 + as) 



{qn - l)Pn - "3 ' 



?34 ■ Kn- Kn = - l)Pn 

+ {ao + as - l)g„ - (02 + a4)(ai + a2 + a4)- 



qnPn - 014 

However 



Tg-/ • - Kn = t{t - 1)- log -^=t{t~ 1)- logr„, 



and we employ the results of l2.30l and the evaluation of the coefficient functions in 
(|3.7I3.8|I to arrive at 



rN UN TN 

{t-lfjL = JjL + N+[N-l + i, + u){l-rNfNf-^. 

rN UN TN 

In addition we note that after recalling (|2.t)|) . H4.11|l factorises into 

{N + 11 + uj){N + ^l + u))rNfN = [qNPN + + i^][{qN - 1)^^ + ^ + ^^]- 
The stated results, (|4.12I4.15|) . then follow. □ 

There is another perspective on r^^ [N] for which the Okamoto r- function theory 
can be used to provide a recurrence system based on dPy distinct from that in 
ProDOsition l4.1l The starting point, used extensively in JI], is to obtain r^^[A^]|^=o 
as specified bv 14.61 via the projection (—1, 1) ^ T of an average over the Jacobi 
unitary ensemble. This relates H4.6|) to the Pyi system with the parameters 

(4.16) (ao, ai, a2, 03, a4) = — ^ — uj, N + 2/i, ~N, —fi — u},N + 2uji^ . 

Sequences of the Hamiltonian variables {qn,Pn, Hn,Tn}n=o,i.,... are now generated 
by the shift operator L^l = r3SiS4S2SoS3S2. It has the action L^l : ai i-^- ai + 
1, a2 1-^ a2 — 1, a4 1-^ a4 + 1. Using the methods of we have the following result. 

Lemma 4.1. The sequence of auxiliary variables {gm fn}n=o.i,... defined by 

(4.17) 5« 

On. - 1 
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(4.18) /„:= 



1 - t 



(qn - t){q„ - l)pn 



+ (1 - ao - a2)(9« - 1) - asiqn ~ t) - 



(qn - t){qn ~ 1) 



generated by the shift operator L-^^ satisfies the discrete Painleve equations 

Proof. Using the action of the fundamental reflections and Dynkin diagram auto- 
morphisms given in Table 1 of we compute the action of L^l on q and write it 
in the following way, 



{q-t){q-t) 



t[q{q — l){q-t)p+{ai+a2)q^-{[ao+ai+a2)t-aa—a4)q-a4t] 



{q-m-l) 

X [q{q - l){q - t)p + [ai + a2)q^ - ((ai + a2)t - ai)q ~ a^t] 
[q{q -l){q- t)p + [ai + a2)q^ - {-a^t + ai + a2)q - a^t] 
[q{q — l){q — t)p + {ai + a2)q^ - (-("s + oi4)t + ai + a2 + as))? - a4t], 

where q :— qn, q '■= q^i+i- From the definitions H4.17I4.18|I this result can be readily 
recast as (jl.fill . The second (|1.7(l follows from a computation for /„ + fn-i using 
the shift operator L14. □ 

Making use of Lemma 14.11 together with elements of the Okamoto theory as 
detailed in ^Jj, and applying this to derive recurrences according to the strategy 
of jni) gives the following recurrence scheme for {t^^[N]}n=o.i,2....- 

Proposition 4.3. Let {gN, fN}N=o,i.---' satisfy the dPy system 

XfN + N + l){fN + N + fi + u;} 



(4.19) 
(4.20) 



QN+igN = t- 



fN + /at-I = M + 



iV + 2/i {N + 2LUi)t 



9n - 1 

where t — e*"^ subject to the initial conditions 

go - t 
go - 1' 



gN - 1 



90 



fo = 



1 



1 - t 



{p. + uj){qn - 1) + (/i + t^)(go - t) - 2lj: 



(go - t){qo - 1) 



go 



with 



(4.21) 



go 



-i—\oge''"^Ti{e"'') 



^ + uj + i—\og e^'"^ri(e^'^) 
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Define {qN,PN}N=o,i,... in terms of {Jn, gN}N=o,i,... by 

(4.22) = 

(4.23) 

ffjv - 1 



Pn 



(1 - 1)9- 



N 



{9N - l)fN + Co)gN + {N + 2lui)'^—^^ ^ N - fi ~ to 

gN ~t 



Then with To{e"f') = 1 and Ti{e"^) = woie"^') as given by {Tn}n=2,3,... 
is specified by the recurrence 

(4.24) -(7V + ^ + ^)(iV + ^ + ^)ZimJ^ 

-'at 

= QNiqN - l)^P% + [(2m - N)qN + N + 2lui] {q^ - 1)pn - ^f^Nq^ - NiN + 2wi). 
Proof. Let y„ := L^^K,, - K„ = K^+i - K„. From Tabic 1 of fTP we have 

{t-l)q„(, IN , , 1 , (1 - ao - a2)(ai + ^2 + 04) 

Yn = — < yq-n - l)p„ + ao + ^2 - 1 + 



g« - 1 I qn{qn - l)Pn + {ctl + a2)qn + 0:4 

Now consider 

^ 1„„ Tn+lT„-l 



tit - 1)- log "+ - Kr^+i + - 2K^, 



This latter difference, upon again consulting Table 1 of turns out to be 
r„ - LuYn = t{t - 1)^ log (^qniqn - i fpl 

+ [(ai + 2a2)qn + cti^qn - l)pn + a2[(Q;i + a2)g„ + 

After integrating both expressions and introducing an integration constant (|4.24(l 
follows. □ 

Remark 4.1. It is known from ^OIiE] that the sequence of auxiliary variables 
{gnj fn}n=o,i,... defined by 

(4.25) 5„ := 

qn -L 

(4.26) /„ := g„(g„ - l)p„ + (1 - ^2 - a4){qn - 1) - aa^n - Qo ^"*''^" ^''^^ 

gri - t 

generated by the shift operator Lq^ satisfy the dPy equations 

t (/n + 1 - a2)(/ri + 1 - ^2 - 04) 



(4.27) gn+ign = 



t-l /„(/„ + as) 

ai agt 



(4.28) /„ + /„-i--a3 - , ■ 

gn - 1 t{gn - 1) - 5„ 
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In fact the two systems of recurrences (14.2714.28(1 and ((1.611.7(1 are related by an 
element of the 5*4 subgroup of the Wa(Fi) transformations, namely the generator 
[221 • This has the action 

(4.29) -.ao^ ai,t^ T~l''^^ T~T''^ '"^ ~ 

and when applying these transformations to ((4.251) . 1(4.26(1 . ((4.27(1 . 1(4.28(1 we recover 

(jIT7|l . (jHHI), iCni), mZIl respectively. 

5. Applications to Physical Models 

5.1. Random Matrix Averages. A specialisation of the above results with great 
interest in the application of random matrices |18 ( .|17 ( is the quantity 

N 

(5.1) F^''^{u■^,)■.= {'Y{\u + z,\^''') . 

This has the interpretation as the average of the 2/i-th power, or equivalently the 
2/i-th moment, of the absolute value of the characteristic polynomial for the CUE. 
In the case |u| = 1 ((5.1() is independent of u and has the well-known (see e.g. [2j) 
gamma function evaluation 

(5.2) lf\\u + z,\^^) -if\\l + zi?') = n 

T=l /CUE„ n=e"* \ ' " /cUE„ f}^ V^{] + l + ll)' 

where for convergence of the integral 5R(/i) > — V^- For \u\ < 1 we see by an 
appropriate change of variables that 

N N 

(5.3) = ^F[^\-^i, -fi;N;ti, . . . , %)|t,=...=t„=|„|2, 

where the second equality follows from 113.50(1 . For \u\ > 1 we can use the simple 
functional equation 

/ = l'"l ' \ I 1 | - + ^'P^/ 

1 = 1 1 = 1 

to relate this case back to the case |m| < 1. 

The weight in the first equality of ((5.3(1 is a special case of ((1.3(1 . In terms of 
the parameters of the form ((1.3(1 we observe that ^ = 0, 2/^ (-^ /i, = (Ii = /i/2, i.e. 
UJ2 ~ and t = The trigonometric moments are 

(5.5) w^n^ iw^^,"^/^ — -2Fi{-^i,-^i + n;n + l;\u\'^) n G Z>o, 

(5.6) Wn = |up"w-„ n e Z>Q. 

The results of Section |31 then allow ((5.1(1 to be computed by a recurrence involving 
the corresponding reflection coefficients. 



(5.4) (ni"+-rL..,H^r-(ni^ 
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Corollary 5.1. The general moments of the characteristic polynomial \ det{u + U)\ 
for arbitrary exponent 2/i with respect to the finite CUE ensemble U G U{N) of rank 
N is given by the system of recurrences 

pCVEpCUE 

(5-7) (•^CUE^2 = 1 ~ l"l ^JVi 

with initial values 

(5.8) F,^^^^l, FP^^ = ,FA-f^,-f^;l;\un 
and 

(5.9) 2\urrNrN-i ~ \u\^ ~ 1 

1 _ |,,|2W 2 

' ' [(^+l + M)|^/pr^+i + (A^-l + M)rAr_i] 



{{N + m)|m|V^v + (iV - 2 + M)riv- 



rN-\ 

subject to the initial values 

5.10) ro = 1, ri = — Tim ■ 

2Fi(-^, -/x; 1; \u\^) 

Proof From either (|^ . (|XTU|) or and the fact that ri = we can 

repeat the arguments of Corollary 13.41 to deduce that fjv = \u\'^^rN for > 0. 
The recurrence relation follows simply from the specialisation of H3.19|l and the 
initial conditions from the = 1 case. □ 

Another spectral statistic of fundamental importance in random matrix theory is 
the generating function E^^{{7t — (p, tt); ^) for the probabilities E^^{k; (vr — 0, tt)) 
that exactly k eigenvalues lie in the interval (tt — 0, tt). This is specified by 

(5.11) ii;cuE((^„^,^);^) 



e ^ — e ''I , 



where Cat = (27r)^A^!, and thus corresponds to the special case fi = uj — id — 0, 
t = e"^ of H1.8|l . We remark that for these parameters the underlying Toeplitz 
matrix elements, given for general parameters by l|3.1() or l|3.2|l . have the elementary 
form 

(5.12) ^i;„=J , 

Zm n 

A recurrence scheme for the generating function E'f^^, deduced as a corollary 
of ProDOsition l4.1l has been presented in |l()j . Here we use recurrences found herein 
for TAT, fN, together with the fact that for l|5.11|l one has rN = t~^fN, to replace 
the role of the coupled recurrences from |1()) by a single recurrence. 
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Corollary 5.2. The generating function for the probability of finding exactly k 
eigenvalues z — e'^ from the ensemble of random N x N unitary matrices within 
the sector of the unit circle 6 G (tt — 0, tt] is given by the following system of 
recurrences in the rank of the ensemble N , 

pCUE pCUE 
W'-L-Jj ^^CUE-)2 ~ 

where the initial values are 

(5.14) E^^^ = l, i?cuE^i_A<^^ 

zvr 

and the auxiliary variables xjq are determined by the quasi-linear third order recur- 
rence relation 

(5.15) 2XNXN-1 " 2 cos ^ = [(^ _^ 1)xn+i + {N - 1)xn-i] 

2 xm 



_ 1 ^jV-l 
XN-l 

or the quadratic second order recurrence relation 



[NXN + {N- 2)xN-2\ , 



(5.16) 

(1 - xlf [{N + lfxl+^ + {N~ lfxl,_^] + 2(iV2 - 1)(1 - x%)xN+iXN-i 
+ 47V cos ^2;jv(l - x^) [{N + 1)xm+i + [N - 1)xn^i] 



2 V 9 

cos' 2 " 



= 0, 



along with the initial values 



(5.17) X-i = 0, xq — 1, xi = 



^1 ^ 



Proof. The first recurrence relation follows directly from the general recurrence 
and CoroUary EH whilst the second follows from (|T77|) . □ 

5.2. 2-D Ising Model. Consider the two-dimensional Ising model with dimcn- 
sionless nearest neighbour couplings equal to Ki and K2 in the x and y directions 
respectively (see e.g. p ). Let ctq.o and crN,N denote the values of the spins at the 
lattice sites (0,0) and {N,N) respectively. For the infinite lattice an unpublished 
result of Onsager (see |23) gives that the diagonal spin-spin correlation (cro,oO'Af,Af) 
has the Toeplitz form 

(5.18) (CTo,oo-Af,Af) = det(ai_j(fc))i<,,j<7v, 

where 



1 



(5.19) ap{k) := — / dOe' 

ZTT 



2tt 



l-(l/fc)e-'^ 
1 - {l/k)e^^ 



1/2 



DISCRETE PAINLEVE EQUATIONS, 



29 



k ~ sinh 2Ki sinh 2/v 2 . The analytic structure is different depending on A: > 1 (low 
temperature phase) or /c < 1 (high temperature phase). 

Jimbo and Miwa JS] identified H5.18|l as the r-function associated with a mon- 
odromy preserving deformation of a linear system, which in turn is associated with 
Pvi . In a more recent work the present authors have identified H5.18|l as a 
r-function in the Okamoto theory of Pyi These identifications have the con- 
sequence of allowing H5.18|l to be characterised in terms of a solution of the u-form 
of the Painlevc VI equation, when regarded as a function of t := k^. 

It is our objective in this subsection to use the recurrences of Subsection 3.1, 
appropriately specialised, to derive a iV- recurrence for H5.18|l . Before doing so, we 
remark that a recent result of Borodin can also be used for the same purpose. 
The latter recurrence applies to all Toeplitz determinants 

(5.20) qi^'^'^?) (1 - 0^^' det[g._,]..,.i,...,„, 
where 

(5.21) g,^^ J^^ dee-^P' (l - ^e^') ' (l - v^e" 

As previously noted in 11 , comparison with 15.18|l and (|5.19|) shows 

oo^ / \ /I ; -2\l/4 (-1/2,1/2,1/fc^) 

(5.22) (cro,oCTw,Jv) = (1 - fc ) Qn ■ 

By regarding the Fourier integral in H5.19|) as a contour integral, and changing 
the contour of integration, the corresponding weight function can be chosen to be 
equal to JI] 

(5.23) zi/4|l + z\-'/^{l + k-^z^/^ = zi/2(i + z)-i/2(l + rt)i/2^ l/fc2 < 1^ 
(5.24) 

fc-lz-3/4|^ ^ ^|-l/2(^ ^ ^2^)1/2 ^ ^-1^-1/2(1 + z)-l/2(i + ^2^)1/2^ ^2 < ^ 

The first, multiplied by /c^/^, is the case 

(5.25) ^ = 0,A* = = c^2 - = l/k^ 
of (|1.3|l . while the second multiplied by fc^/^ is the case 

(5.26) ^^0,/i-i,^i = -J,W2 = -^,i = A:2. 

Substituting these values in (|3.1|l gives Gauss hypergeometric evaluations for the 
matrix elements (|5.19|) . After an appropriate limiting procedure to account for 
factors which otherwise would be singular, one obtains the well known fact that the 
Toeplitz elements in the low temperature regime are given by 

(5.27) ^ ^^""^ '^f\['^'\ F,{-y2,n+i/2;n + l;k-'), n > 0, 

TT 1 (n -f- 1) 

/ -iyi+lt^-2ri p/ _ l/2lrf3/2') 

(5.28) Wr,= ^-^ iin /^M Pi ^j.^(y,^^„y,.„^i.fc-2)^ „ > 0^ 

TT 1 (n + 1) 
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whilst those in the high temperature regime are 

,5.2<„ . V-':';:p' .i^.(V.~-fA:n + 2;.»). „ > 0. 

TT i (n + 2) 

(5.30) ^ ''i'^^^'^'^ 2fi(-V2, n - 1/2; n; fc^), n > 0. 

The resuhs of Subsection 3.1 provide the following recurrence scheme. 

Corollary 5.3. The diagonal correlation function for the Ising model valid in both 
the low and high temperature phases (with fc i— s- 1/fc m the latter case) is determined 
by 

/ro1^ (o'0,00'A'+l.JV+l}(o'o,oO'JV-l,Ar-l) , 

(5.31) -5 = l-rArrjv, 

along with the quasi-linear 2/1 

(5.32) {2N + 3)k-^{l - rjvrjv)?-iv+i + IN [fc^^ + l - (27V - l)fc-Vjvrjv-i] rjv 

+ (2^ - 3) [(2iV - l)rjvfjv + 1] ^at-i = 0, 

and 1/2 recurrence relations 

(5.33) (2A^ + 1)(1 - rNfN)fN+i + 2iV [(2iV - 3)fArrAr_i + fc'^ + 1] fj^ 

+ (2iV - l)fc-2 [_(2Ar + l)rjvfAr + 1] fAr_i = 0, 

subject to initial conditions for the low temperature regime 

fc2-lK(fc-i) k^-lK(k-'^) 

(5.34) .0 = fo = l,n = ^ + -^^^ = -1 + ^ eH' 

or to the initial conditions for the high temperature regime given by 
(5.35) 

1 r 2 E(k) ] fc^Effc) 

J-o ^-0 = 1, ri = - <—r - , , , > , ri = 



3 ifc2 (fc2 _ i)K(fc) +E(fc) J ' ' (fc2 - l)K(fc) +E(fc)' 
where K(fc),E(fc) are f/ie complete elliptic integrals of the first and second kind 
respectively. 

Proof. The equations (|5. 3215. 33(1 follow from H3.33|l and its "conjugate" upon spe- 
cialisation of the parameters as required by H5.25|l and (|5.26|) . The initial conditions 
follow from (|3.21l) . (|5.27|) - (|5.3U|) and the relationship between the Gauss hypergeo- 
metric functions therein and the complete elliptic integrals. □ 

A consequence of the results of Subsection 3.3 is the following result relating 
the quantities of Corollary 15.31 to the generalised hypergeometric function 2^1^^- 

Corollary 5.4. In the low temperature phase the diagonal correlation function is 
given by 

(5.36) (cro,oCrjV,Jv) =2-P'l^^(-V2,V2;^;il,---,*w)|ii = ...=t« = l/fc2> 
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whilst the reflection coefficients are given by 
(5.37) 



rN 



(5.38) 



^! A^\-y2,y2-N-M,--.,tN) 



^{N - l)\ lim,^o 62-^^1 H-V2, - 1 + e; ti, 



=tjv = l/fc2 



(V2) 



JV 



2Fi(^'(-V2,V2;A^;ii,...,iM) 



tl = ... = tjV = l//c^ 



In the high temperature phase the diagonal correlation function is 
(5.39) (ao,o^Ar,Ar) = I^^5^^fc^2^(i)(i/2,i/2;iV + l;ti,...,i^)|t^. 
and the reflection coefficients are given by 

,,,„, , ,,A- ;fi"'('A,V2;A' + 2;ti,....liv) 



-tj\T — k^ ■) 



(5.41) 



\A 



, in) 



(V2)iV 2^(1) (1/2, 1/2; iV+l;ii,...,tAr) 



l = ... = tjV = fe2 



Proof. The evaluations (|5.36|l and (|5.39(l follow by specialising the parameters in 
(|3.50|l as required by H5.25|l and (|5.26(l respectively. The evaluations 1)5.37(1 , H5.38|l in 
the low temperature phase follow from (|3.51|) , (|3.52|) by specialising the parameters 
as required by H5.25() . Some care needs to be taken with fjv because —fi + uj = 0. 
Applying a limiting process leads to 



(5.42) lim e^F^^^ (-1/2, - 1/2; - 1 + e; ti 



,tN) 



E 



tN 



k:1{k)=N 



(1) 



[N-iy. 



in which only those terms with lengths 1{k) ~ N contribute to the sum. The 
high temperature expressions follow from the low temperature ones through the 
transformation /i «-> o^i . □ 

It is of interest to note that as N grows more of the leading order terms in 
the expansion of (|5.36|) become independent of N, and the following limit becomes 
explicit 



lim {ao.o<JN.N) - (1 - k-Y^\ 

N~^oo 



(5.43) 

At zero temperature, k — 00, the solutions simplify to 

Tat = (iV > 1), (cro,oCTA,A) = 1, 



(5.44) r, = (-l)-(-'/^)- 
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whilst at the critical point, k = 1, we have the simple solutions 



(5-45) rN= fN = {-ir, In 

N 



(2iV + 1)(27V- 1)' ' ' ' 2N + V 

N 

(5.46) {o-o^oo-n,n) = Y\. 



^.^^ro- + V2)r(j-V2)' 

At infinite temperature they become 

(5.47) = (-1)^1^^, = (-1)^^, (ao,oa^,Ar) = (iV > 1). 
These are all in agreement with the known results (21, . 
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